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We consider the stochastic dynamics of Ising ferromagnets (either pure or random) near zero 
temperature. The master equation satisfying detailed balance can be mapped onto a quantum 
Hamiltonian which has an exact zero-energy ground state representing the thermal equilibrium. 
The largest relaxation time teq governing the convergence towards this Boltzmann equilibrium in 
finite-size systems is determined by the lowest non- vanishing eigenvalue Ei = 1 /teq of the quantum 
Hamiltonian H. We introduce and study a real-space renormalization procedure for the quantum 
Ifamiltonian associated to the single-spin-flip dynamics of Ising ferromagnets near zero temperature. 
We solve explicitly the renormalization flow for two cases, (i) For the one-dimensional random 
ferromagnetic chain with free boundary conditions, the largest relaxation time t^q can be expressed 
in terms of the set of random couplings for various choices of the dynamical transition rates. The 
validity of these RG results in d = 1 is checked by comparison with another approach, (ii) For the 
pure Ising model on a Cayley tree of branching ratio K, we compute the exponential growth of 
teq{N) with the number A*' of generations. 

I. INTRODUCTION 

The stochastic dynamics of classical Ising ferromagnets has been much studied for fifty years [H, In particular, 
many works have been devoted to the domain growth dynamics at low temperature T < Tc (or at zero temperature 
T = in ff = 1 where the critical temperature vanishes Tc = 0) when the initial condition is random (see the review 
on phase ordering dynamics [5]). 

In the present paper, we do not consider the dynamics starting from a random initial condition, but we focus instead 
on the largest relaxation time t^q needed for a finite systems to converge towards thermal equilibrium. This largest 
relaxation time teq is defined as the inverse of the smallest non- vanishing eigenvalue Ei of the time-evolution operator. 
Near zero-temperature, more precisely when the temperature is much smaller than any ferromagnetic coupling Jij 

< T < (1) 

the thermal equilibrium is dominated by the two ferromagnetic ground states where all spins take the same value, 
and the largest relaxation time t^q corresponds to the time needed to go from one ground state (where all spins take 
the value +1) to the opposite ground state (where all spins take the value —1). We should stress that we consider 
that the temperature is arbitrarily small, but does not vanish, so that the transition between the two ground states 
is possible and the final state of the dynamics is unique (For studies on the zero-temperature dynamics, where the 
spin-flips corresponding to an energy-increase become impossible, we refer to the recent works [J| and to references 
therein) . 

Of course near zero temperature, the equilibration time t^q becomes extremely large, and numerical simulations of 
the microscopic dynamics become inefficient. Here we introduce and study a real-space renormalization procedure 
valid near zero temperature for the dynamics of pure or random ferromagnets. The paper is organized as follows. In 
section ini we recall the standard mapping between the master equation describing the stochastic dynamics of classical 
systems at temperature T and a special type of quantum Hamiltonians that have an exact zero-energy eigenstate. In 
section IIIIl we describe the various choices of dynamical transition rates for single-spin- flip dynamics of the classical 
Ising model, and the corresponding quantum Hamiltonians. In scction lTVl wc introduce the real-space renormalization 
procedure for these general quantum Hamiltonians. In section [Vj we show how a closed RG procedure can be defined 
and exactly solved for the random ferromagnetic chain. In section IVTl we solve the RG fiow for the pure Ising model 
on a Cayley tree. Section [VIII summarizes our conclusions. In Appendix |X1 we describe another approach that allows 
to check the validity of the RG results of section [V] Finally in Appendix [B] we discuss the contributions that can 
depend on the choice of transition rates satisfying detailed balance. 
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II. RELAXATION OF CLASSICAL MODELS TOWARDS THERMAL EQUILIBRIUM 
A. Master Equation satisfying detailed balance 

To define the stochastic dynamics of a classical system, it is convenient to consider the master equation 

dPt (C) 



Y,Pt{C')W{C' ^C)~Pt{C)Wout{C) (2) 
c 

that describes the time evolution of the probability Pt{C) to be in configuration C at time t. The notation W {C C) 
represents the transition rate per unit time from configuration C to C, and 

Wont {C) = Y,W{C^ C) (3) 
c 

represents the total exit rate out of configuration C. 

For a classical system where each configuration C has some energy U{C), the convergence towards Boltzmann 
equilibrium at temperature T = in any finite system 

Pe,{C) = (4) 

where Z is the partition function 

Z^^e-^^^'^) (5) 
c 

can be ensured by imposing the detailed balance property 

g-/3t/(c)^ (C ^ C') = e-^^^^'^W (C C) (6) 



B. Mapping onto a Schrodinger equation in configuration space 

As is well known (see for instance the textbooks 0-01 )j the non-symmetric operator describing the stochastic 
dynamics of a classical model at temperature T can be transformed into a symmetric quantum Hamiltonian problem. 
In the field of disordered systems, this mapping has been much used for one-dimensional continuous models [8Hll[, 
and more recently for many-body spin systems like the Sherrington-Kirkpatrik model ([ 12] and Appendix B of (131). 
In the field of pure spin models, this mapping has been used for more than fifty years 

In the present context, this standard mapping consists in the change of variable 

P,{C) ^ e~i^^''>MC) = e-^^^'^' < C\yJt > (7) 
Then the master equation of Eq. [2] becomes the imaginary-time Schrodinger equation for the ket > 

^IV* -H\i;t > (8) 

where the quantum Hamiltonian 

n = ^€{C)\C ><c\+J2 ^(<^' >< C\ (9) 



contains the on-site energies 



e(C) = I^™t(C) = ;^I^(C^C') (10) 



and the hoppings (using EqlB]) 

Vice') - -e-t[^(c')-c/(c)]^(C' ^ c) = -e-^[^('^)-^('^')]iy (C ^ C) 

= -V'wlc'TcyW^^lC^^rC) (11) 
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C. Properties of the quantum Hamiltonian "H 

Let us note En the eigenvalues of the quantum Hamihonian H and |i/)„ > the associated normalized eigenvectors 

n\^n > - > (12) 

c 

The evolution operator e~*^ can be expanded on the eigenstates 

e-t^^Y^e'^'^'li^nXi^nl (14) 

n 

The conditional probability (C|Co) to be in configuration C at i if one starts from the configuration Cq at time 
t — can be written as 

Pt {C\Co) - e"f [^(c)-c/(Co)] < c|e-*«|Co >= e-^[^('^)-^('^")l ^ e-^"V«(C)V:(Co) (15) 

n 

The quantum Hamiltonian H has special properties that come from its relation to the dynamical master equation : 

(i) the ground state energy is Eq = 0, and the corresponding eigenvector is given by 

lV'o>=E^^|C> (16) 

The normalization l/\/Z comes from the quantum normalization of Eq. [T3l 

This property ensures the convergence towards the Boltzmann equilibrium in Eq. [7] for any initial condition Co 

P, (C|Co) ^ ^ e-^[^('^)-^('^«)]Vo(C)V^(Co) = = Fe,(C) (17) 

(ii) the other energies i?„ > determine the relaxation towards equilibrium. In particular, the lowest non-vanishing 
energy Ex determines the largest relaxation time {1/Ei) of the system 

Pt{C\Ca)-Peq{C) ~ e-^i*e-^[^(^)-^(C")lVi(C)Vr(Co) (18) 

Since this largest relaxation time represents the 'equilibrium time', i.e. the characteristic time needed to converge 
towards equilibrium, we will use the following notation 



tea = ;^ (19) 



In summary, the relaxation time t^q can be computed without simulating the dynamics by any method able to 
compute the first excited energy Ex of the quantum Hamiltonian % (where the ground state is given by Eq. [12] and 
has for eigenvalue Eq = 0). For instance in [l^], the 'conjugate gradient' method has been used to study numerically 
the statistics of the largest relaxation time in various disordered models. Let us now describe more precisely how this 
general framework applies to single-spin-flip dynamics of Ising models. 

III. QUANTUM HAMILTONIAN ASSOCIATED TO SINGLE-SPIN-FLIP DYNAMICS 

A. Single-spin-flip dynamics 

We consider a system of classical spins Si = ±1 where each configuration C = Si, S2, ■■ has for energy 

U{C) = -Y.J.jS,S, (20) 

i<j 

The couplings Jij may be random. 
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Within a single-spin flip dynamics, the configuration \C >= \Si > {82 > --ISn > containing A'' spins is connected 
via the transition rates W {C ^ C) to the TV configurations obtained by the fiip of a single spin Sk —Sk denoted 

by 



\Ck >=a^k\C> 

in terms of the Pauli matrix 172:. The energy difference between the two configurations reads 

U{Ck)-U{C)=2SkY,Jk^S, 

i 

The quantum Hamiltonian of Eqs [^ [TUlfTT] can be thus rewritten as 

N 

n =Y,eiC)\C><C\ + J2J2 ^(C, a^,C)at\C >< C\ 



N 



C fe=l 
!l[UiC^)-U{C)] 



,-§lU(Ck)-U{C)] 



\c ><c\ 



(21) 



(22) 



(23) 



C fc=l 



B. Simplest choice of the transition rates 

It is clear from Eq. [231 that the simplest quantum Hamiltonian corresponds to the following choice of the dynamical 
transition rate 



W{C = {S,} ^ Cu = alO = e-4[£/(c.)-c/(C)] ^ g-/?5.(E.^. ^..^0 
The quantum Hamiltonian of Eq. [23] then reads in terms of Pauli matrices 



(24) 



v. 



simple 



C fe=l 
N 

EE[^ 



-|[C/(Cfc)-C/(C)] _ 



C k=l 
N 



|C ><C| 
\C ><C\ 



fe=i 



(25) 



where we have used the identity 1 = \C >< C|. The quantum Hamiltonian of Eq. 1251 has been mentioned as the 
simplest for the one-dimensional ferromagnetic chain in Eq (4) of Ref [isj . 



The Glauber choice for the transition rates Q 
WiC^ {S^} ^ Ck = <jIC) = 



Glauber choice 



^-!^[U(Ck)-U(C)] 



2 cosh (I [U{Ck) - U{C)]^ 2 cosh (/3 [^.^fe J^kS^ 
corresponds to the more complicated quantum Hamiltonian 

1 



H 



Glauber 



N 

E 



1 2 cosh 



/3 {J2^^k ■hk^Tf 



(26) 



(27) 



where we have used the fact that cr^ has for eigenvalues (±1) and that cosh is an even function. 

This quantum Hamiltonian of Eq. [571 has been used already used for the Sherrington-Kirkpatrick spin-glass model 
and for the finite dimensional ferromagnetic Ising model in Ref. [l3| (see Appendix B and Appendix C respectively). 
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For the one-dimensional pure ferromagnetic chain, where each spin Sk has only two neighbors, the local field Bk = 
JikSi = J{Sk-i + Sk+i) can take only the three values h^. = — 2 J, 0, 2 J, so that one may replace the exponential 
factors using projection operators to recover the forms given in Refs [15l - [l7l | . This type of 'first quantized' quantum 
spin Hamiltonian can be transformed further into 'second quantized' Hamiltonian involving annihilation/creation 
operators or Fermi operators using Jordan- Wigner transformation (see the review (isj ) : this method has been 
followed in particular by Ref. for the Glauber dynamics of the pure Ising chain. However in the present paper, 
we will work directly on the 'first quantized' form of the quantum spin Hamiltonian of Eq[2S]or Eg 1271 with the aim 
to define a real-space renormalization approach, in analogy with the Strong Disorder RG rules introduced for the 
random transverse field Ising model on its first-quantized form (see the discussion in section flV Al below) . 

D. Most general choice 

To better understand the structure of the renormalized Hamiltonian that will be generated by the real space RG 
procedure introduced in the following section, it is useful to consider the most general choice satisfying the detailed 
balance equation of Eq. |6] 

WiC = {S^} -^Ck = <jIC) = Gk{Su .., Sk-i,Sk+u SN)e-"^''^'^^*^ ^^"^^^ (28) 

where Gk{Si, .., Sk-i, Sk+i, ••, Sn) is an arbitrary positive function of the (N — 1) spins i ^ k that may depend on 
the index k. The positivity requirement 

Gk{Si, .., Sk-i, Sk+i, ■-, Sn) > (29) 

ensures that all elementary single-fiip are possible with a non- vanishing rate, so that the ground state of Eq. [16] is 
unique, and the dynamics converges towards thermal equilibrium. Since the reversed transition of Eq. [55] has the 
following rate 

W{Ck ^alC^C = {S,}) = Gk{Si,..,Sk-uSk+u...SN)e+^^'^\^^*^'^^^'\ (30) 

the amplitude Gk{Si, Sk~i, Sk+i, Sn) represents the 'symmetric part' of the two opposite transitions involving 
the fiip of the spin k (when all other spins remain fixed) 

Gfc(5i,..,5fc_i,5fe+i,..,57v) = {C ulC)W {alC ^ C) (31) 
The corresponding quantum Hamiltonian of Eg reads 

N 

^general = ^ Gfc (cT^ . . , i , a^^l , • • , ) (e"^-^ (^M. '^--r ) - a.^) (32) 
fc=l 

In finite dimensions with only nearest-neighbors interactions, it seems natural to consider local rates where the 
amplitude only involves the spins i that are neighbors of k (i.e. such that Jki 7^ 0) 

G'fc(crJ, ..,crfc_l,crfc+l, ..,cr^) ^ G''k""''' {Wi}jik^o}) (33) 

To respect the symmetries of the classical energy of Eq. 1201 a further requirement is usually that the amplitude should 
only be an even function G{x) = G{—x) of the local field 

In the next section, we introduce a renormalization procedure for the amplitude Gk and we will see that even if one 
starts with a function of the local field (Eq. [34)) . this form is not stable via renormalization and will generally lead to 
a local function of the neighboring spins (Eq. [55]) . 

IV. RENORMALIZATION RULES NEAR ZERO TEMPERATURE 

In this section, we derive the appropriate renormalization rules for the quantum Hamiltonian introduced in the 
previous section in relation with the single-spin fiip dynamics of classical ferromagnets. 
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A. Differences with the disordered quantum Ising model 

Let us first mention that in the high-temperature limit /3Jij = ^ 1, the quantum Hamiltonian of Eq. [55] or [27] 
reduces to the standard transverse-field Ising model 

jjsrraple,Glauber ^^^^^^ ^ ^ ^ ^^f^^. j^^^.^ _ (35) 

When the coupling Jij are random, the low-energy physics is then well described by the Strong Disorder RG procedure 
valid both in one dimension [19] and in higher dimensions d > 1 [2Qi-|2^ (see [23j for a review). However here we are 
interested into the opposite limit of very low temperature where /3Jij — '-^ ^ 1 (Eq[T]), where one cannot linearize 
the exponentials in the quantum Hamiltonians of Eq. [55] and 1271 In the following, we thus derive appropriate RG 
rules in this opposite regime. 

B. Analysis of an elementary operator 

The general quantum Hamiltonian of Eq [31] can be considered as a sum 

N 

^general ^ ^ (3g) 
fc=l 

of elementary operators 

^ Gk{al, .., a^+i, .., a^) (e-^<^^ ) - at) (37) 
The corresponding matrix elements in the cr^ basis 

< S[, .., S'MSi,.., Sn >= '5s^.,5, j Gfc(5i, .., Sk-i.Sk+u Sn) < Slle^^'^'^i^^" '^-^0 _ at\Sk > (38) 
are diagonal for all spins j ^ k. 

1. Effective problem for a single spin 
For each fixed value of the local field Bk = X)j JikSi, we may diagonalize the effective problem for the single-spin k 

hlff{B,)^e-f'<^''-al (39) 
The two normalized eigenstates are respectively 

\vk{Bu)> ^ 1^ J e^^-|^fc = +l>+e-|^'=|^fc = -l>j (40) 

A/2cosh(/3Bfc) L J 



and 



with eigenvalues 



\wk{Bu)> ^ \ e-^">'\Sk = +l>-e^">'\Su^-l> (41) 

^2cosh(/3Bfc) 



hl^^{Bk)\vk{Bk)> =0 

hl^f{Bk)\wk{Bk)> =2cos\i{PBk)\wk{Bk)> (42) 
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so that the single-spin hamiltonian of Eq. |39]can be rewritten as the projector 

'■k 



hl"(.Bk) = 2cosh{(3Bk)\wk{Bk) >< Wk{Bk) 



e-^^'l^fe - +1 > -e-^'=|^fc = -1 >J ^e-^^" < ^fc = +1| - e^""" < Sk = -1| 

= ale-^^^^'' [\Sk - +1 > +\Sk = -1 >] [< Sk = +1|+ <Sk^ -1|] a|e-4^^^'= (43) 
In terms of operators, the elementary operator of Eq. 1371 may be thus rewritten as 

hk = Gfc(cri, .., ct|_i, cr|+l, ..,cr^) 

a'^e-^'^i J.,<^} = +1 > +1^^, = _i >] [< Sk = +1|+ < Sk = a^^e"^-^ (44) 

^. Properties of elementary operators 

Eq|33]is convenient to see explicitly the positivity property for any ket \ip > 
<ip\hk\ip> = ^ GkiSi,..,Sk-i,Sk+i,:,SN) 

5i,..,5fc-l,5fc + i,..,5Ar 

|e-4 ^'■.^.^(^^^ _^ 5fe_i, 5fc - +1, 5fe+i, .., 5w) - ^- -^-^'^(^i, .., 5fe = -1, Sk+i, .., 5^)P 
> (45) 

as a consequence of the positivity of the Gk (Eg [^^ . Moreover, it is clear that the exactly known ground state of 
zero energy of Eq. [16] 

|V-o> = <Si,..SN\iJo>\Si,..,SM>= J—- e^^i<.<.<"^-^'^^-|^i,..,5jv > 

Zn{/3) = ^ e/3Ei<.<j<iv''.jS,s, (-^g^ 

Si,..Sn 

satisfies 

'0o(5'i, .., Sk — +1, .^fc+i, .., Sn) _ ^0j2i JkiS, ^^y^ 
tpoiSi,.., Sk-i: Sk = —1, Sk+i, .., S'at) 
and is thus annihilated by all elementary operators /i^ for fc = 1, .., 

hk\^o>^0 (48) 

as it should. 

C. Renormalization of the sum of two neighboring elementary operators 

1. Sum of two neighboring elementary operators 

The sum of two neighboring local operators of Eq. [37] of index k and I with Jki ^ reads 
hk + hi =G,(aJ,..,a^_i,a^+i,..,a^)(e-^-^(^--°+^.^'^--')-afc-) 

+G,(a?,..,af_i,af+i,..,a^) (e-/'-?(^--^+E,v.^'.-j) -af) (49) 
The corresponding matrix elements in the cr^ basis 

<Sl..,S'j^\hk + hi\Si,..,SN> =( n ^s'^'S, 

[Gk{Si,..,Sk-i,Sk+i,..,SN) < ^^|e~^'^^(''-^'+^'^'''-^0 > 
+Gi{Si,..,Si-i,Si+i,..,Sn) < 5;'|e-'^'^'(^-^^+^^^'=-'-^') -afl^z >] (50) 

are diagonal for all spins j ^ (fc, I). 



8 



2. Effective two-spin problem 
For each fixed value of all the other external spins Sj^(k,i), using the notations 

Bk = ^ Jki Si 



Bi ^^Jo^, (51) 

and 

fff' = Gfc(5i, .., Sk-i: Sk+1, •, 5*;., Sn) 

gf"" =Gi{Si,.,Sk,--,Si-i,Si+i,..,SN) (52) 
we have to diagonalize the effective problem for the two spins {k,l) 

= gl' (e-^^^^-^--' - at) + gf' (e-^-r(^--i+sO _ af) (53) 



The four-dimensional vector 



l'^A>= E cf'^'|5fe,5,> (54) 

Sfe = ±,S,=± 



is an eigenvector of the operator h'^f/ of Eq. [53] with eigenvalue A if 

= [e-^'^- {gte-^^^ + g^e'^^^) X] C++ - 3+0+- - .9+<:,-+ 

= [e^-^^' {g^e-^^^ + g^e^^') X] 4" - ,+4+ - 5.7c,- 

= [e'^-^- (g+e^^'= + 5,-^-^^') - A] c,-+ - .g+c++ - g^c^ (55) 

5. Finding the lowest non-zero eigenvalue Ai 
We already know the exact ground state |ua=o > of eigenvalue A = 0, with components (not normalized here) 

To find the exact three other eigenvalues, one has to solve the remaining cubic equation. Here, since we are interested 
in the other small eigenvalue Ai <^ 6^'''=', we will neglect Ai in the two last equations of Eq. [55] to obtain 

r+- = e-^J"' 3t4^+9kcr 

^ g+ePB.+g-e~PB, ^^'^ 

that we may replace in the two first equations of Eq. [S5] 

U - Lff, e +g,e X,e \ c, g, ^ ^^^^^^^ g, ^^^^^^ ^ 

-k7e/^^+gre^"'-A,e/^^-]c,--.g,7 '^Ss^f^^.B. - 9^ l'Sj'l%s. (58) 

gj, e P"'' + g^eP^' g^eP"'' + gi e p"' 

For Ai — 0, we recover the exact solution of Eq. [56] as it should. The other eigenvalue reads 

Ai - e "/!costi /5(i:ffc 3—— — , _ — + + arR -i-R -I , — , „ ^ — (59J 

9t9k +9U1 +9U^ef'^'''+^'^ +9k9i e-^i^.+B,) 



9 



with the corresponding components of the eigenvector \uxi > (not normalized here) 

(•++ — g§Jki-§Bk~§Bi 

The two other components c^^ and c'^^ are given by Eg 1571 

4- Projection onto the two ferromagnetic states 
The projection of the operator of Eq. [S3] onto its two lowest states of eigenvalues Aq = et Ai reads 

K7 ^ ^ ^ \ux, XuaJ (61) 



At the level of approximation we are working, we wish to keep only the two ferromagnetic states ++ and (the 

two other states H — and — h have been taken into account in Eq. [57] to produce renormalized rates between H — h and 
in Eg I58p . So we keep only the two following leading components in the eigenvector 

\ux, >~ e^'"' (e-^(^^-+^')| + + > -e^iB,+BO\ _ _ >^ (62) 

with the corresponding normalization 

< ux, \ux, >^ e^'^'='2 cosh[/3(Bfe + Bi)] (63) 

Eq [HT] becomes 



2cosh[/3(Bfc + B0] 

^ ^ > _ _ _ >^ |^^-|(B,+B,) < ^ ^ I _ JiB,+B,) < _ _ |) (64) 

So the two spins Sk and Si now form a single renormalized ferromagnetic cluster, that we may represent by a single 
spin with the two states 

\+>R = I + + > 

\->R =\--> (65) 

with the renormalized external field (see Eg [5T|) 

Br = Bk+Bi = (Jkr + Ju)S^ (66) 

that corresponds to the natural renormalization of ferromagnetic coupling between the new cluster {kl) and i 

JOcl)^ = Jk^ + JU (67) 

In terms of this renormalized cluster, Eq. 1641 reads 

hl^/ ^ Gn (e-^^«|+ >n -e^^-\- >n) (e-^^- < +\h - e^^« < -|«) (68) 

It has the same form as an elementary operator of Eq. [31] with the renormalized amplitude (using Eq. 

Ai 



2cosh[l3{Bk + Bi)] 

„-l3Jki ^fc ^fe y-n - • ^ I ' ill m xnk^ • nk^ i 

^ +-,+-, 4- + «R„ , - - _«R„ V"^-* 



at 9k (.9,+e^^^ + .9i e-^ ^'^)+g+g, (ff,+e^^' + e '^^Q 



9k 9k +9i 9i + 9k 9i e'^"« + 9k 9i e 
in terms of the variables defined in Eq. [S2]for a given value of the external spins Si^(^k,i) 
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5. Renormalization rules for operators 

At the level of operators in the whole Hilbert space, the renormalized Hamiltonian of Eq. |64] reads (using the same 
notations as in Eq. 137] ) 

h^i) ^ G^/K, af+i., a^) (e-^-«^« - aj,) (70) 

where the amplitude of Eq. |69] is the renormalized operator 

G..(a„..,a,_„a,^„.a,_„a,+,,a,) -e g^g- + g^g^ + gtg^e^^^'^^B, + g-g-,-,,B.^B, (^l) 

in terms of the operators 

9k = G'fc(crJ,..,cr|_i,cr|+i,.,S'; =±,.,cr|^) 
5; = G;((Ti,.,S'fe ±,..,CTf_i,crf+i,..,cr|^) 

6. Final State of the RG procedure 

For a system of N spins, we will obtain after {N — 1) RG steps a single renormalized spin Sr — ±1 representing 
the two ground states of the whole sample, where all spins take the same value Sr. Since there is no renormalized 
external field left Br = 0, the final effective Hamiltonian for the single spin Sr simply reads 

^fu^al ^ Qfn^al (^y) (1 - ffjj) (73) 

where G-'^™''' ( A^) is a numerical amplitude. The quantum ground state of zero energy is the symmetric combination 
of the two classical ferromagnetic ground states as it should 

l^^ = + >^^^ = -> (74) 

whereas the excited quantum eigenstate is the antisymmetric combination of the two classical ferromagnetic ground 
states 

= + > -1^^^ - - > (75) 
v2 

with the eigenvalue 

eI"""' (N) ^ 20^"""' (N) (76) 
The conclusion is that the equilibrium time of Eq. \W\ can be thus obtained from the final renormalization amplitude 

Lt^ do 

= ;^7w^ = 2Gf'^^'{N) ^^^^ 

In summary, we have thus defined a renormalized procedure for the amplitudes G that allows to obtain in the end 
the equilibration time teq{N). Now to better understand the meaning of the RG rules of Eq. [TU let us give some 
explicit examples. 
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7. Example : first renormalization step starting from the simple Hamiltonian of Eq. 

The quantum Hamihonian of Eq. [25] corresponds to the initial simple case where all amplitudes are unity 

GT^^h ■•7CTfc„i,CTfc+i, --jcr^) = 1 (78) 

The first RG step where k and / are grouped into a single renormalized cluster yields the following amplitude fEg 1711 
with = 1^ gf 

qR _ ^-PJu, cosh(/3Bfc) + cosh(/3B0 



^kl 



l + cosh(/3(Bfe+Bi)) 



l + cosh(/3E^^(,0{^fe. + ^H)^f) 

So this renormalized amplitude does not remain a number as in the initial condition of Eg [751 but becomes an operator 
that involves the neighboring spins of k and I. 

8. Example : first renormalization step starting from the Glauber Hamiltonian of Eq. \27\ 
The quantum Hamiltonian of Eq. [57] corresponds to the initial case 

GrK,..,aLi,a^+i,..,a^) = 2cosh(/3E. ^^"^ 



so that the operators of Eq. [72] reads 



± 1 



2 cosh (3{Bk±Jki) 

9t =:T— TIJtWtTT (81) 



2cosh/3(B,± JfeO 
Eq[7T]for the renormalized amplitude then becomes 



'kl 



2cosh/3(Bfc + Bi) + 2e2-''«'+'B^--Bi 2cosh/3(Bfe + Bi) + 2e2-''fei-'Bfc+s, 

1 

2cosm3Y.^^^k,l){Jk^ + + 2e2'^-+^'^('=.o(''--^-)-?) 

1 



(82) 



2cosh(/3E,_,(fc ;)(Jfe, + JhX) + 2e2'^'''-^.^(-.o('^'='-'^'')'^f) 

So this renormalized amplitude does not remain of the Glauber form of Eq. [80] in terms of the renormalized local 
field Br — Bk + Bi. Actually it does not even remain a function of the single renormalized local field Br. 

9. Discussion 

These two examples show that the renormalized amplitudes Gr become generically a function of the cr^ operators 
of the neighboring spins, ie of the form of Eq. [33] It is not clear to us at this stage how to determine the operator 
form that would remain stable upon the general RG rule of Eq. 1711 In the following, we will thus concentrate on two 
geometries where we can obtain closed RG rules, namely the one-dimensional case and the Cayley tree. 

V. BOUNDARY RENORMALIZATION FOR THE RANDOM FERROMAGNETIC CHAIN 

In this section, we consider the random ferromagnetic chain of N spins with the classical energy (Eq I20p 

U{Si, ...Sn) = ~ JiA+iSiSi+i (83) 



with free boundary conditions for the two boundary spins Si and Spf. The couplings J^^^+i are positive random 
variables. 
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A. Closed RG for the simple Hamiltonian of Eg 1251 in one dimension 

1. First step of the RG rules in the bulk 

The first RG step where the two neighboring sites (fc, fc + 1) are grouped into a single ferromagnetic cluster leads 
to the following renormalized amplitude (Eq[7Hl) ^oi 1 < k < N ~ 1 

qR ^ ^-/3Jfe fe^i C0sh(/3Jfc_i,fc) + COSh(^Jfc+i,fc+2) 

Since there remains the renormalized local field operator [ Jfe_i.fe(j^_j^ + Jk+i,k+2<^k+2\ denominator, this renor- 

malized amplitude is not just a number as the initial condition. To avoid this difficulty, we will now consider what 
happens near one boundary. 



2. First step of the RG rules near the bow 

We now consider the case fc = 1, where we make a ferromagnetic cluster out of the two sites (1,2) near the boundary. 
Since the spin zero does not exist, we have Jk-i,k — Jo,i = , so that Eq. 1841 reduces to a renormalized number 
(without any operator anymore) 

Gf, ^,-/^.,. l + cosh(^J2,3) ^^-,.,,, 
' 1 + cosh(/3J2,3) 

This shows that we may define a simple closed boundary RG procedure as follows. 



3. Closed Boundary RG procedure 



If we iterate the renormalization near the boundary, the renormalized state after (n — 1) RG steps will be the 
following : the spins (l,2,..,n) have been grouped together into a single renormalized spin with some associated 
renormalized amplitude ^j. The other spins {n + 1, N) of the chain are still in their initial form with amplitude 
unity Gfc = 1 for n + 1 < fc < A'^. Let us now perform the n RG step where the spin (n + 1) is included in the boundary 
cluster. The RG rule of Eq. [71] yields the recurrence 



riR 



^+1] 



_„+ir-<R 



2G(^„] +2cosh(/3J„+i 
iGfl,n])^ + 1 + G«„](2cosh(/3J„+i,„+2))_ 



(86) 



Within the bulk, i.e. for > 0, in order to be consistent with our previous approximations in the low- 

temperature limit (Eq[T]), we may replace 



2cosh(/3J„+i_„+2)) ^ ^ 6' 



/3J„ + l,„ + 2 



(87) 



Moreover, we expect that the renormalized amplitudes G^ only decay upon renormalization, and are thus smaller 
than their initial unity value (see already the first step of Eq. I55|) 



G, 



l,n] 



< 1 



So Eq. 15^ reduces to 



riR 



-/3J„,„ + i/^-R 



p/3^n + l,ji + 2 



(89) 



i.e. this recurrence becomes simpler in terms of inverse variables 
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4- Final result for a chain of N spins 
For a finite chain of N spins where Jn,n+i — 0, the recurrence of Eq. |86] yields for the last step with n + 1 — N 



Rfinal _ -/3Jjv-l w/^^J 

^[l,Af-l] 



n-- — p 



2^[l,JV-l] + 2 



-PJn-i,n r^R 



'-'[l.Ar-ll 



2 



1 + ^[l,Af-ll 



i.e. using inverse variables 



1 1 



^Rfinal o 
[1,JV] 



□/3 Ja 



r<R 

^[l,Af-l] 



We now now use iteratively the recurrence of Eq. |90] valid in the bulk to obtain 



1 



r<R 

^[l,N-2\ 



g/3JAf-l,JV _j_ ^0JN~2,N~l + f^JN-l,N _j_ g/3JjV-3,iV-2 + /3JjV-2,iV-l 



o/3 Jjv-3,JV-2 

"-"[l.TV-S] 



1 ^ 



:,^C'^fc-l.fc + <^/c,fc+l) 



(91) 



(92) 



(93) 



(with the notations Jq.i = = Jn,n+i) 

The conclusion of this RG procedure is thus that the equilibration time of the finite random chain of N spins reads 
for the 'simple' dynamics (Eq. [77|) 



1 



N 



(94) 



B. Closed RG for dynamics depending only on the local Held in one dimension 

We now consider the more general case where the amplitude Gk are a single even function G{x) — G{~x) of the 
local field (see Eq. [M]). 

(95) 



Gk - G (J/c-i,feCTfe_i + Jk.k+i^k+i) 
The Glauber Hamiltonian of Eq [57] corresponds to the special function 

G{x) ^ 



2 cosh(^a;) 



(96) 



1. Closed Boundary RG procedure 



If we iterate the renormalization near the boundary as in section IV A 31 we obtain finally the following recurrence 
(Eqini) : 



in terms of the numbers 



g[^,n](/»+l + + /++l/-+l(2C0Sh(/3J„+i,„+2)) 

(G« „,)2 + + G« „](/++ie'5''"+i."+^ + /-^^e-'5^"+i,"+^) 

fn+l = G {Jn. n+l ^ J n+l, 11+2) 



(97) 
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Within the bulk, i.e. for Jn+i,n+2 > 0, in order to be consistent with our previous approximations in the low- 
temperature hmit (Eq[T]), we may use Eqs 1871 and the fact that the renormahzed amphtudes only decay upon 
renormalization to simplify Eq [97] into 



~ „-/3J„,„+i^iJ 

^[l,n+l] - e Lz[i, 



Jn+lJn+l'-- 



i.e. this recurrence becomes simpler in terms of inverse variables 

g/3Jn + l,ji + 2 g/3J„,„ + l+/3J„ + l,„ + 2 gfiJn,n+l 



r<R 

^[l,n+l] 



fn+1 



(99) 



(100) 



2. Fmal result for a chain of N spins 
For a finite chain of N spins where Jn,n+i — 0, Eg l98l reads 

Jn = G ( Jat-i^at ± 0) = /jv 

and the recurrence of Eq. [97] yields for the last step with n + 1 ^ N 



" '-'fi.w-il 



^—0Jn-i.n/~<R 

^[l,Af-l] 



Vn 



N 



i.e. using inverse variables 



1 



1 



y^Rfinal o 



^PJn-i.n g/3Jjv-i,jv 



[1,^-1] 

We now now use iteratively the recurrence of Eq. llOOl valid in the bulk to obtain 



r<R 

[1,JV] 



r<R 



JN 

^ 0{Jk-i.k + Jk.k + i) 



iN-2 



k=l 



fk 



r<R 

""[l,Af-3] 



(101) 



(102) 



(103) 



(104) 



(with the notations Jo,i = = Jn,n+i) 

Using Eq. [35] the conclusion of this RG procedure is thus that the equilibration time of the finite random chain of 
N spins reads for the dynamics defined by the amplitudes of Eq. l97l 



1 



^ g/3(./fc-i.fc + ./fc,fc + i) 



i?/"-'(7V) 2G^_^] 4 ^ G {.h-i,k - JkM+i) 
In particular for the Glauber dynamics corresponding to Eq. I96[ the equilibration time reads 

N N 
tglauber^j^^ ^ g^(.7.-l.. + + l ) 3 cosh (/3 Jfc_i,fc - /3Jfe,fe+i) - - ^ [g2/J J.-i.. + g2/3J.,. + l] 



(105) 



fc=i 



fc=i 



N-l 



1 + E 



(106) 



fe=l 
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C. Discussion 



In this section, we have obtained a closed RG procedure in dimension d — 1 that yields the explicit result of Eq. 
llOSI for the equilibration time of a finite chain. The validity of this RG result is checked in Appendix |^ using another 
method. Since we compare the explicit expression in terms of the set of all random couplings, for various choices of 
the transition rates, this agreement shows the exactness of the RG procedure near zero temperature. 

The physical meaning of Eq. 11051 is that the equilibration time is the sum of N random variables, possibly sligthly 
correlated (the same random coupling Jn.n+i appears in the two terms corresponding to fc = n and k = n + 1). For 
the Glauber case of Eq. 11061 even these slight correlations disappear and the equilibration time reduces to the sum 
of independent random variables e^^'^'-' ^'+i whose distribution can be computed from the distribution of the random 
couplings Jfc,fc+i. 

Note that in the limit of a pure chain where all ferromagnetic coupling have the same value J, Eq ll05l reduces to 



G(J) ' 'G(O) 



« -gW 



as expected : the Arrhenius factor e^'^'' comes from the barrier (2 J) to create a domain-wall at one boundary, and 
the prefactor TV comes from the small probability of order 1 /N that this domain- wall crosses the whole system of size 
N instead of returning back. For the Glauber case where G(0) = 2, Eq. 11071 is in agreement with the leading term 
near zero temperature of the open pure chain discussed in Refs [13] ■ 



VI. BOUNDARY RENORMALIZATION FOR THE ISING MODEL ON THE CAYLEY TREE 



In this section, we consider the pure ferromagnetic Ising model with the classical energy of Eq[20] defined on Cayley 
tree of branching ratio K (coordinence (K + 1)) with N generations, and with free boundary conditions on all the 
boundary spins. We focus on the dynamics corresponding to the simplest choice of Eq. 1241 where the corresponding 
quantum Hamiltonian of Eq. [25] has initial amplitude Gk all equal to unity. 

We wish to define a closed boundary RG procedure that preserves the symmetry between the K offsprings of a given 
site. So the basic RG step concerns K renormalized boundary spins (S'l, 5*2, .., 5*;^) whose renormalized dynamics is 
described by some renormalized amplitude G (which is a number and not an operator) and their common ancestor 
spin S whose dynamics is still described by the initial amplitude 5 = 1 and by the external field Ba = Ja^a induced 
by its next ancestor spin Sa ■ The ferromagnetic couplings have still their initial value J so that we have to study the 
following effective Hamiltonian for these {K + 1) spins (^i, .., Sk, S) 



jjsirnple _ |^g-;3.^ (Ef., J-!+B.) _ j + g (^-^J-' _ (108) 

i=l 

After the first RG step, the renormalized amplitude G is expected to become smaller and smaller, so the most 
appropriate approach is a perturbative analysis in the parameter G. 



A. Properties of the Hamiltonian Hj^^^ for G = 

For G = 0, the spins (Si, .-iSk) cannot fiip and are thus frozen, so the problem reduces to the single spin S in 
the external field B = {^2!i=\ '-^^i + ^^^^ have already studied in section HVB II So the 2^+^ states can be 
decomposed into 

(i) the 2^ states corresponding to Eq. HHl 

|^5„S.,..,5. ^ ^J]!^. >E / ' l'^> (109) 

J = l S=±^ 2 cosh PiY^f^^ J S^+Ba) 

that have a vanishing eigenvalue for G = 
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The physical interpretation is that the spin S is at equilibrium with respect to the frozen spins {Si, .., Sk)- 
(ii) the 2^ states corresponding to Eq. Hi] 



j=i s=± V2cosh/3(^^i JS, + Ba) 



that have a finite eigenvalue for G = 



K 



2cosh;3(^ JS, + Ba) 



w, 



Note that the exact ground state fEq \W\i of the Hamiltonian of Eq. 
1 



IV'o > = 



V Zk+1 



Si=±S2=± Sk=±S=± 



^K,. ^ E E E E - 

Si = ±S2 = ± Sfc=±S=± 



Si=±S2=± Sfc=± 

belongs to the subspace spanned by the 2^ states \vq^'^^'"'^ 



^'"^^ > > ...\sk >\s> 

K 

2 cosh JS^+Ba) 

Sfc=± L i=i 

> fEalTTUI 



\^o> =^E E- E A 



A' 



2cosh/3(E + Sa)|wo^ 



> 



(111) 



(112) 



(113) 



(114) 



B. Perturbation in the parameter G 

We have seen that for G — 0, there are 2^ states (Eq lllOp that have zero energy. For small G > 0, the perturbation 
will lift this degeneracy : the exact ground state of Eq. 11141 will keep its zero energy for arbitrary G, but the 
other (2^ — 1) eigenvalues will become positive as soon as G > 0. To determine them, we need to diagonalize the 
perturbation within the subspace spanned by the 2^ vectors - Sa- • g_ look for an eigenstatc via the 

linear combination 



\ux> ^ E ^ 

Si=±,S2 = ±,..,Sk=± 



Si,S2,---,Sk\,Si,S2,---Sk 

ro 



> 



(115) 



The eigenvalue-equation 



=(i/^';^f'^-A)|uA> 



^ ^ 2^5'i,S'2,---,5/i- 



EG(e-^^?^^^-<)-A 



Si ,S2 5 ■ ■ ■ 5ft 



> 



(116) 



can be projected onto the 2^ bra < y^^'^^'-'^'i^ ^ to obtain a system of 2^ linear equations for the 2^ coefiicients 

yS'i,S2,---,Sjr 



=<vf-''-\iH'ji^f''-X)\ux> 



=< V, 

From the matrix elements 

S -I So T - S 



S[,..S'i( I 



E 



S2,---,Sk 



EG(e-''-?^-^-ar)-A 



Si,...Sk 



> 



(117) 



5i,S2,---5ft: 




>= ri'^^^-. 



2cosh/3(^^._,.^ J^, +B,) 



2cosh/3(^;ti JS'j+S,) 

-m^.s 1^. s 2cosh/3(E,,.J^.+i?.) 

J " V2cosh/3(E,^, - JS, + Sj^2cosh/3(E,^, ^S", + JS^ + B,) 



(118) 
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we obtain that Eg 11171 yields 



= 

Si = ±,...,Sk=± 
K 



,Sk 



K 



^2cosh/3ef^i JS'^+Ba) 



/jnS;„...,Sj^ 

R ^ 

- T'^i'--'~'^^'---''^g 



2cosh/3(Ej^, JS'^ - JSi + Ba) j2cosh/3(E,^, JS'^ + JS[ + S,,) 



(119) 



Let us now use the symmetry between the K branches to note t{k) the components T^^'^^'-'-'^ii where k spins take 
the value (+) and {K — k) take the value — , i.e. 



(120) 



For the two extremal cases t{K) (ah spins are (+)) and t{0) (all spins are (— )), Eg 11191 becomes 

2cosh/3((/<i:- 1)J + BJ 



= 
= 



• 2cosh/3((X-l)J + i3a) _ 

2cosh/3(i^J + BQ) 

• 2cosh/3((JC-l)J-i3a) _ ■ 

2cosh^(i^J-B<,) 



t{K) - GK 
t{Q) - GK 



V2 cosh P{{K ~2)J + Ba)^2 cosh l3{K J + Ba 
2cosh^((is:- l)J-Ba) 



-.t{K-l) 



y^2coshl3{{K - 2)J - Ba)^/2coshl3{KJ - Ba) 



t{l) (121) 



whereas for the non-extremal case < k < K , Eq. 11191 reads 
= 



^^^ 2cosh/?((2fc^ j^- l)J + 5a) ^ _ ^^^ 2cosh/3((2fc-i^ + l)J + g,) _ ^ 



2cosh^((2fc-ii:)J + Ba) 



-Gk- 



2 cosh P{{2k ~ K -1)J + Ba) 



^2cosh/3((2fc - K - 2)J + B„)^2cosh/3((2fc - K)J + Ba) 
2 cosh /3((2fc - K + 1)J + Ba) 



-G{K~k) 



2 cosh P {{2k - K)J + B a) 
t{k-l) 

t{k + l) 



m 



a/2 cosh /3((2fc -K + i)J + Ba)y/2coshl3{{2k - K)J + Ba) 
It is easy to check that the components (not normalized here) of the exact ground state of Eq 11131 



io(fc) - ^J2cosh|3{{2k - K)J + Ba) 



(122) 



(123) 



satisfy Eouations 11211 and 11221 for A = as it should. 

To compute the lowest non- vanishing eigenvalue Ai, it is consistent to set Ai = in all egs 11221 concerning the 
non-extremal cases {) < k < leading to 



= 



,^ 2cosh/?((2fc-A'-l) + 5a) ^ _ ^^ 2 cosh/3((2fc - K + I) + Ba) 

y/2 cosh I3{{2k - K) J + Ba) 
2 cosh P{{2k - K - 1) + Ba) 



^2coshl3{{2k - K)J + Ba) 



m 



^2 cosh/3((2fc -K -2) + Ba) 



t{k -1)-{K- -) 2cosh/3((2fc-i. + l) + i^.) ^ „ ^ 
\/2cosh/3((2fc-ii:4-2) + Ba) 



that may be solved in terms of the boundary conditions t{K) and t{0). 

To obtain the explicit solution, it is convenient to introduce the amplitudes A{k) with respect to the ground state 
components of Eq. [13 



t{k) = A{k)to{k) = A(fc)v/2cosh/3((2fc - K)J + Ba) 
so that Eq. 11241 takes the simpler form 

A{k) = p^{k)A{k - 1) + p+{k)A{k + 1) 



(125) 
(126) 
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with the notations 

fc2 cosh^((2fc - K - 1) + Ba) 



P-{k) 



[/fc2cosh/3((2fc - K - 1) + Ba) + (K ~ fc)2cosh^((2fc - K + 1) + Ba)] 



^rj,. = (i^-fc)2cosh/3((2fc-j^ + l) + i3„) = 1 - Wfc) (127) 

' ~ [k2coshp{{2k- K -1) + Ba) + {K-k)2 cosh I3{{2k-K + 1) + Ba)\ ^ ' 

Let us introduce two hnearly independent solutions. The sohition corresponding to the boundary conditions 

QxiO) =0 

Qk{K) - 1 (128) 
^ RKil,k) 



can be obtained by recurrence [25| and reads 



using Kesten variables [26| 

i?K(l,0) -0 
Rk{1,1) =1 

k—lrn / N 

RK{i.k>2) =i+^n^^ 

^-^ v+(n) 



R a K) 1 I V TT^-^") 1 1 P-^^^ I P-Wp-(^) p_(l)p_(2)...;._(i^-l) 
Similarly, the solution corresponding to the boundary conditions 

Qo(0) - 1 

Qo{K) =0 (131) 

R,{k,K-l) 
^i?o(0,i^-l) 

Ra{K,K-l) =0 (133) 
i?o(-f^- 1) =1 

i?o(fc<i^-2,i^-l) =1+ Jj'^+W 



reads 



with 



v^in) 

m—k+\ n=m 

p+{n) ^ , P+{K-l) , , p+{K-l)p+{K-2)...p+{l) 



P-{n) P-(J^ — 1 

t-n — 1 n — ->-)■) ^ \ / ^ \ / 



=i„=iP-(") P-iK-l) p-(i^-l)p-(if-2)...p_(l) 

It is useful to introduce the continuation of the ground state components of Eq 11231 to half-integers to rewrite the 
ratios 

p_(fc) _ A:2cosh/3((2fc-i^-l) + S,) _ kt^ {k - ^) ^^^^^ 



p+{k) {K-k)2coshp{{2k-K+l) + Ba) {K -k)tl{k+\) 



and the products 



m / X m 



ntl {n-\ 



2) 



m\{K-l-m)\ tl{\) 



_{K-n)tl{n+\)\ {K-l)\ (m + i) 



m 2^ 
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In particular in the following, we will need the two denominators 



A'-l 



Rk{1,K) 



il{K-l-m)\ tl{^) 



7n—l 



(K-iy. (g(m + i) 



1 



K-1 

2V ^0 C^-itU^+l) 



and 



RoiO,K-l) 



{m - ly.jK - my. tl {K - \) 



that determine the solutions near the boundaries of Qk 

Rk{1,1) _ 1 



= ti{K- 



K-l 

E 



(136) 



(137) 



and of Qo 



Qk{1) 
1 - Qk{K - 1) 

Qo(i^-l) 
i-Qo(i) 



Rk{1,K) Rk{I..K) 
^ Rk{1,K)-Rk{1,K-1) 
Rk{1,K) 

Ra{K -l,K _ 1 



i?o(0,A'- 1) 



(138) 



RoiO,K-l) RoiO,K-l) 
Ro{0,K-l)-Ra{lJ<-l) _ 1 



Ro{0,K-l) Rk{I.K) 
The solution of the system [T26l that satisfy the boundary conditions fEg I125P 

m 



(139) 



A(0) = 
A{K) = 



io(0) 
t{K) 
to{K) 



can be obtained by the linear combination 



A{k) = A{0)Qo{k) + A{K)QK{k) - ^Mk) + PjlT^QKik) 



toiO) 



toiK) 



so that the solution of the system 11241 reads 

t{k) ^ A{k)to{k) ^ 



Qa{k) + -m7\QK[k) 



To determine Ai, we just need to replace 

"^o(l) 



t{l) = 
t{K-\) = 



Lto(0)^°^^) 



m 



toil) 



to(fc) 



(140) 



(141) 



(142) 



toiO) 



QoiK-1) 



toiK) 
tiO) 



'.K 



(1) 



tiK) 



tojK-l) 
toiK) 



QxiK-l) 



tiK) 



(143) 



in Eas ll21l to obtain the following system of two linear equations for the two components iAi(O) and tx-^iK) 

tliK) 



= 







1-Qa-(/^-1)-Ai 



GKtl {K^\) 



txAK)-*j^QoiK-l)t,M 



1-Qo(1)-Ai 



tm 



txAO)~^^QKii)txAK) 



GKtl {I) 

The two components (not normalized here) are orthogonal to Eq 11231 as it should and read 

tx, (0) = toiK) = ^2 cosh PiK. J + Ba) 
tx, iK) = -to(0) - -v/2cosh/3(ifJ-B,) 



(144) 



(145) 
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The corresponding eigenvalue Ai reads using Eqs I136[ I137[ [T38l Eg 11391 



GK 



^K-l 1 



C™__j2cosh,9((2m+l-_fS')J+B„) 



2 cosh P{K J - Ba) 2 cosh I3{K J + Ba) 



(146) 



Since this expression is unchanged via the transformation Ba — JSa — > —Ba, we may replace Ba by its absolute value 
\Ba\ = J to obtain the final expression for the lowest non- vanishing eigenvalue Ai at first order in perturbation with 
respect to the parameter G 



Ai 



GK 



m=0 



C"^_-^2 cosh 13 J(2m+2-K) 



1 



2 cosh p J {K -I) 2cosh/3J(i4: + l) 



0{G^) 



(147) 



C. Renormalization rule for the amplitude 

Let us now project the Hamiltonian of Eq. 11081 onto its two lowest eigenvalues Aq = and Ai 



H 



simple 
K+1 



Ai 



|uAi >< UaJ 



(148) 



where the eigenvector of Eq. 11151 can be approximated at low temperature by its two components onto fully ferro- 
magnetic states 



K 



K 



y/2 cosh I3{K J + Ba) 



\s> 



\s > 



y/2 cosh 13{K J - Ba) 

tx, (K) I n = 1 > 1 1^ = +1 > +tx, (0) I n = -1 > 1 1^ = -1 > 



K 



(149) 



with the coefficients fEq ll45p 



tAi(O) = ^ 2 cosh P{K J + Ba) 



\(^K.J+Ba) 



Finally at leading order near zero temperature, one obtains 

K 



y/2cosh(3{KJ- Ba) ~ -e2(-^'^--s» 



KJ 



n = 1 > 1 15 = +1 > -c 



with the corresponding normalization (using Ba — JSa) 



< UaJuAi > 

In terms of the renormalized spin 



fS{KJ+B^) _^ ^I3{KJ-B^) _ J3{KJ) 



Y[\s, = -\>\\s = -i> 



e'"^^"^2cosh(/3Ba) = e'^(^^)2 cosh(/3J) 



(150) 



(151) 



(152) 



K 



\Sh = -> ^ (ni^, = -l>| |^ = -> 



(153) 
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and of the external local field Bu = Ba = JSa, the effective Hamiltonian of Eq. [148] can be rewritten as an elementary 
operator (Eq. 



H 



simple 
K+1 



■){' 



e^^« < Sf. 



with the renormalized amplitude (using Eg I147P 

Ai GK 



Gi 



2 cosh /3 J 2 cosh S 7 V''^"^ r 



2cosh^J(is: - 1) 2cosh/3J(ii: + l) 



(154) 



(155) 



To be consistent with the previous low-temperature approximations, we now should evaluate the leading behavior 
of Eq. 11551 near zero temperature, i.e. we should replace hyperbolic functions by exponentials. In particular, one has 



1 



1 



1 -f ^1 

2 '^ = Y-^ 



2 cosh pj{2m + 2~K) e/3"'(2™+2-A:) + g-/3 J(2m+2-K) 
so that the leading term near low temperature of Eq. 11551 depends on the parity of K 



K 



^ 2 



(156) 



1. Leading behavior near zero temperature for even K 

When the branching ratio K is even, then (f - 1) is an integer, so that the integer m can take this value, and the 
sum in the denominator of Eq. 11551 is dominated by this contribution 

K-l ^ ^ 

2 C|?_i2cosh/3J(2m + 2-if) ^ ^^^^^ 



so that Eq 11551 reads at leading order 



K even : G« ^ Ge-^'''2KCj:l = Ge-P-'^2 ^ ^ ^ (158) 



For instance for K ~ 2, one obtains 

K = 2 : G" ~4e~2'^"'G (159) 



2. Leading behavior near zero temperature for odd K 

When the branching ratio K is odd, then (-^ — l) is not an integer so that m cannot take this value, and the sum 
in the denominator is dominated by the contributions of the two closest integers to = ^Z^^ and ni 



K~l 
2 



K-l 

E 



^-f}.J\2m+2-K\ 



^0 C^'_i2 cosh/3 J(2to + 2-K) 



E 

m=0 



-I3J 



c 



K-l 



C 2 C 2 



{K-iy. [K-iy. 

-,j K{^y{^y 



so that Eq 11551 reads at leading order 



K odd : G^ ^ Ge-^'^^^-^^ ^ Ger^'^(^'^\K-l)cft2 



For instance for X = 3, one obtains 



(160) 



(161) 



K 



(162) 
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D. Conclusion for the equilibrium time t''^^^^''{N) of a Cayley tree with N generations 

Let us now consider a finite Cayley tree of brancfiing ratio K witli N generations. For tlic first RG step where 
Go — 1, we cannot use the perturbative analysis presented above to obtain Gi. However since Gi ^ 1 at low 
temperature, we may use the perturbative analysis given above to obtain the recursion 

^ ^ ^ j^hr ^^^^^ 

for all RG steps corresponding 1 < n < — 1, where the factor p{K) has been evaluated at low temperature (Eqs 
[TSglandfTMD 

p^K) = ^ ^ = e^^^Ui_iA^ for even K 



Finally, at the last RG step, we could take into account that the center has {K + 1) neighbors instead of AT, and has 
no further ancestor Ba = 0. However, this anomalous last step is only a boundary multiplicative contribution, as is 
Gi, and cannot change the dependence upon the number N of generations for large N coming from Eg 11631 

1 

[p(k)] 



^n""" « 771^ (165) 



In summary, we obtain that the equilibrium time ig*™P''^(iV) = l/(2Gj^*"°') (Eg [771) of ^ Cayley tree of branching 
ratio K grows exponentially with the number N of generations 

i-™P'^(7V)cx (166) 

where the growth factor p{K) is given explicitly by Eg 1 1641 for any K. 



E. Equilibrium time t^,'''"'"=''(iV) for the Glauber dynamics 

The above results concerning the simple dynamics can be extended to the Glauber dynamics as follows. The 
Hamiltonian H^^f''^ of Eq. 11081 has to be replaced for the first step by 

2cosHp{j:'Lj-!+Ba))^ ^ 2cosh(/3Ja^)^l ^) ' 

Since the K leaves have no external field and are just linked to a^, the amplitude 2 coshog jo-^) reduces to the number 
2cosh(^j) ■ "^^^ remaining non-trivial amplitude - — h(ff(X]f 1 Ja-'+B )) ^^^^ disappear when we apply the perturbation 
method within the subspace annihilating the corresponding operator (^e~^°' (5Di=i +b°) _ u^^. Our conclusion is 

thus that the equilibrium time t^'""'"^'" (TV) for the Glauber dynamics will have exactly the same leading exponential 
behavior in N as the result of Eq. 11661 derived for the simple dynamics 

even if the prefactor can differ (see the discussion on the differences between the equilibrium times of the simple and 
Glauber dynamics in Appendix |B]) . 

F. Comparison with previous results on dynamical barriers 

From Eq. 11661 and Eq. 11681 we obtain that the energetic barrier Bk{N) defined as the coefficient of /? in \nteq{N) 

Bk{N)^ lim "'^ o = NJK + 0{1) for even K 

= NJ{K-1) + 0{1) for odd if (169) 
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grows linearly with the number A'' of generations (i.e. logarithmically with the number of sites Mn oc K^) in 
agreement with previous works of physicists (27l - l29| and of mathematicians [30l - [33j . Besides this correct scaling with 
N, it appears that the slope (A'— 1) J for odd K of Eq. II 691 coincides with the slope obtained in [i^, where a so-called 
'disjoint strategy' is optimal, whereas the slope KJ for even K of Eq. 11691 differs from the slope J{K — 1) obtained 
in [23, H^, where a so-called 'non-disjoint strategy' is optimal. We refer to Refs [l^-dH for more explanations on 
the differences between disjoint/non-disjoint strategies. For the present work, it is clear that the renormalization 
procedure making coherent clusters of spins within sub-trees corresponds to the disjoint strategy. 

Finally, besides the Arrhenius factor involving the energetic barrier of Eq |1691 the present renormalization procedure 
predicts explicit combinatorial prefactors for the exponential growth factor p(K) (Eq. I164p that have not been 
previously discussed in the literature, to the best of our knowledge. 



VII. CONCLUSION 



In this paper, we have introduced a real-space RG procedure valid near zero-temperature to evaluate the largest 
relaxation time of classical random ferromagnets. We have used the standard mapping between the master equa- 
tions satisfying detailed balanced and quantum Hamiltonians having an exact zero-energy ground state. The largest 
relaxation time teq governing the convergence of the dynamics towards the Boltzmann equilibrium is determined by 
the lowest non- vanishing eigenvalue Ei = 1/teq of the quantum Hamiltonian H. We have thus defined appropriate 
real-space RG rules for the quantum Hamiltonian to evaluate Ei for finite systems. We have described how the 
renormalization flow can be explicitly solved for the two following cases. 

(i) For the one-dimensional random ferromagnetic chain with free boundary conditions, the largest relaxation time 
teq can be expressed in terms of the set of random couplings for various choices of the dynamical transition rates. 
The validity of these RG results in d = 1 have been checked by comparison with another approach in Appendix. 

(ii) For the pure Ising model on a Cayley tree of branching ratio K (coordinence {K + 1)), we have computed the 
exponential growth of tea(N) with the number N of generations. 

In a companion paper |34l |. we explain how the renormalization flow can be also explicitly solved for the Dyson 
hierarchical Dyson Ising model. In the future, we hope to obtain numerical results for the RG flow in finite dimensions 
d>l. 



Acknowledgments 

It is a pleasure to thank G. Semerjian for his comments on our work and for pointing out the references (29l-[35 



Appendix A: Check of the validity of the RG procedure in d = 1 



In this Appendix, we present another approach to check the results of the RG procedure obtained in section |V] for 
the random ferromagnetic chain 



1. Ansatz for the first excited quantum state in terms of exit probabilities 



a. Eigenequation for ipi 

For the quantum Hamiltonian H corresponding to the spin-fiip dynamics of classical spin models with the energy 
of Eq. [201 the exact ground state of Eq. [TBI 

MC) - ^ (Al) 



is invariant under a global flip of all the spins 

M-c) = MC) (A2) 
On the contrary, the first excited state will be antisymmetric under a global flip of all the spins 

M-c) = -MC) (A3) 
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but its modulus is expected to coincide nearly with ipo{C) in the two valleys around the two classical ferromagnetic 
ground states. It is thus convenient to set 

MC)^MC)A{C) (A4) 
and to look for the antisymmetric amplitude A{C) (antisymmetric under a global flip of all the spins) 

A{^C) ^ A{C) (A5) 
The eigenvalue equation for the quantum Hamiltonian of Eq. 191101 1111 

HlV^i >=^i|V^i > (A6) 

becomes via the change of variables of Eq. IA4I 

[WoutiC) - El] A{C) = Y,W{C^ C')A{C') (A7) 

c 

For a large system where Ei is small, we expect that Ei can be neglected with respect to Wout{C) for all configu- 
rations different from the two classical ground states, so that one obtains the approximate equation 

where 

_ WiC~,C') WjC^C) 

represents the probability that the first exit from configuration C leads to C" for the master equation of Eq. [2l with 
the normalization 



J2nc{C') = l (AlO) 



6. Relation with exit probabilities 



Exit probabilities are known to satisfy backward master equation similar to Eq. IA8I (see for instance the textbooks 
P-Q)- More precisely, in a ferromagnet, one may introduce the probability Q+{C) that the dynamics starting in 
configuration C reaches first the configuration C+ (all spins plus) than the configuration C_ (all spins minus). The 
complementary probability Q-{C) = 1 — Q+(C) represents the probability that the dynamics starting in configuration 
C reaches first the configuration C_ than the configuration C+ . The escape probability satisfies the backward master 
equation 

Q+{C) = Y,^c{C')Q+{C') (All) 

c 

for all configurations C different from the two ground states, and the boundary conditions 

Q+{c+) = 1 

Q+(C_) =0 (A12) 
This suggests the following Ansatz for the antisymmetric A{C) satisfying Eq. lASI up to a normalization factor M 

j^ansatz^(j^ = M {2Q +{C) - 1) = U {I - 2Q_{C)) (Al3) 

using Q+{C) = O-(-C) one obtains A(-C) = -A(C). 

The only point where Q+{C) does not satisfy Eq. lAllI are the two boundaries C+ and C_ where Qj^ is given by 
the b.c. corresponding to 

^ =2Q+(C+)-l = l 

j^ansatz ^ 



AT 



2Q+(C_)-1 = -1 (A14) 
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Let us now estimate Ei for the Ansatz of Eq. IA13I corresponding to 

^Ansatz^(j^ = V'O (C) A'^""'*" (C) = A/'?/'^""''*" (C) (2Q+ (C) - 1) 



(A15) 



via 



El = 



^ ^Ansatz ^^^^^Ansatz ^ 
^ ^Ansatz Ansatz \, 



(A16) 



In the numerator, ah configurations C different from C+ and C_ give zero-contributions as a consequence of Eq. lAllI 
So the only contributions in the numerator come from C = C+ and from C = C_ where we may use the boundary 
conditions of Eq. IA14I to obtain 



E. 



Wout{C+) - 5]VK(C+ ^ C'){2Q+{C') - 1) 



Wout{C-) - yW{C- ^ C")(l - 2Q+(C')) 



c 



EcV'o'(C)(2Q+(C)-l)2 



(A17) 



The numerator involves the probability to reach first C- before returning to C+ when one leaves C+, and the 
probability to reach first C_ before returning to C+ when one leaves C+ , which are the same by symmetry. 



2. Application to the random ferromagnetic chain near zero temperature 



We now focus on the random ferromagnetic chain of N spins of Eq with free boundary conditions for the two 
boundary spins 5*1 and Sm- Near zero temperature (Eq[T]), we may neglect the configurations containing more than 
one domain-wall, and work within the space of the following (2A^) configurations 

\k >'j^"' - [\Si = .. = Sk = -1; Sk+i = .. = Sn = +1> +\Si = .. = Sk = 1; Sk+i = .. = Sn = -I >] 

\k >^^^" = [\Si = .. = Sk = -l;Sk+i = .. = Sn = +1> -\Si = .. = Sk = l:Sk+i = .. = Sn = -I >\A18) 

where A: = 0, 1, ...,N — 1. In physical terms, |0 >''i'™ and |0 >''^w'" are the symmetric and antisymmetric combination 
of the two ferromagnetic ground states where all spin have the same signs, whereas \ksym > and \kasym > with 
1 < A: < A'^ — 1 are the symmetric and antisymmetric combination of the states where there exists a single domain- wall 
between the sites (fc, k + 1). 

We consider the quantum Hamiltonian 



n 



N 



'l^'^l{jk-l,kTl-i+Jk,k+lcrl + i) fff 



2<k<N-l 

where the amplitudes Gk are given by a single even function G{x) — G{—x) of the local field (see Eq. 



G (jN-lM<^k-ll 



(A19) 



a. Two first eigenvectors within the single domain-wall approximation 
The ground state \tpo > of zero energy is exactly known from Eq. lAll 



Sn > 



Si,...,Sj\ 



(A20) 
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Within the reduced space of configurations containing no more than one domain- waU (Eg IA18|) . the ground state 
reduces to 



N-l 



fe=l 



near zero temperature 

To respect the antisymmetry of Eq. IA3[ the first excited state wih be a hnear combination of the antisymmetric 
states of Eq. [XTSl 



>N~ |0 >%'y"' + E e~^J''-^+'AN{k)\k (A22) 



k=l 

with some amplitudes Aj\j(k) that we wish to determine. The eigenvalue equation for this first excited state of the 
Hamiltonian of Eq. IA19I reads 

0- (Hat-^i)!^! >Ar (A23) 

= |0 >^^^" [-E, + f-e-^'^'^ (1 - A^il)) + fj,e-^''--^-- (f + A^iN - 1))] 
+ |1 >^^^™ [-E.e-^'^'^A^il) + /f {An{1) - 1) + f^e-^'^^^ (Ajv(l) - Ar,{2))] 

+ E Ifc >r'" [-E^e-^^'-^^^ANik) + f^e-f"'"-^^" {Ai,{k) - A^ik - 1)) + f^+,e-^'^^^^>'+- (A^ik) - Ai,{k + 1))] 

k=2 

+ \N - 1 >^^^" [-i?ie-'^-^«-^'"^jv(iV - f ) + /^_ie-^^"-'"-^ {An{N - 1) - An{N - 2)) + {An{N - 1) + f )] 
in terms of the numbers 

= G {Jk-i.k ~ JKk+i) (A24) 

b. Ansatz with exit probabilities 

Instead of solving exactly the eigenvalue problem of an iV x iV matrix of Eq. IA231 we have proposed in section 
I A li the following approximation : in all coefficients involving \k >^^™ with k = \, ..,N — \, we may neglect the term 
containing Ei with respect to the others to obtain the {N — 1) equations for A: = 1, .., — 1 

y-e-/3J.-l,. I^^ansatz^^^ _ ^ansatz _ -^^^ ^ j-^^g-/3 J,+i,,+2 (^™«at^(y!j) _ + 1)) = (A25) 

with the following boundary conditions 

y4^""°*^(7V) = -1 (A26) 

The only remaining term in Eg IA23l is then the first line involving |0 >^^'" that determines the value of the energy 
El as 

^ansatz^^^ = f^e-^^^-''- (A^"""*^ (0) - A^''"^' [l)) + fj^e-^^" {Af^''''' {N - 1) - A^^'^'^'iN)) (A27) 

So we have replaced the eigenvalue problem of Eq. IA23l bv a simpler homogeneous recurrence equation (Eg |A25]) with 
the boundary equations of Eq. IA261 that can be solved as follows. 

c. Exact solution for exit probabilities in one dimension 
It is convenient to set as in Eq. IA13I 

= 2Qo(fc) -1 = 1- 2QNik) (A28) 

where QAr(fc) satisfies 

QN{k) ^ P+{k)QN{k + 1) + p-{k)QN{k ~ 1) (A29) 
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with the respective probabilities 

P+{k) 
P-(k) 

and the boundary conditions 



QNiO) =0 

Qn{N) = 1 



(A30) 



(A31) 



Then QN{k) represents the probabihty to reach first the boundary k = N rather than the boundary k — for a 
random walker starting at k and moving with probabilities of Eq. IA30I The well-known solution of this standard 
problem can be obtained by recurrence [23 using Kesten variables ^2^] and reads 



(k) 



R{l,k) 



(A32) 



with 



i?(l,0) 

R{l,k> 2) 
i?(l,7V) 



k—1 rn 



-"--^ P+{n) 

i — l n— 1 ' ^ ^ 
N-1 m , ^ 

V-(n) 



, P-(l)p-(2) 



p_(l)p_(2)...p_(jV-l) 
p+(l)p+(2)...p+(A^-l) 



The corresponding estimate of the energy of Eq. IA27I reads using Eq. IA28I 

^a„.at.(^) = /f e-'^^i.^ ((1 - 2QAr(0)) - (1 - 2QAr(l))) + /^e"^-^"-!^" (1 - 2Qn{N - 1) - (1 - 2Qn{N))) 
= 2f-e-^J''- (QAr(l) - QAr(O)) + 2/^e-'5-^«-i'« (gAr(A^) ~QNiN~ 1)) 
_ ^/re-^'^v2i?(i, 1) + /^g-/J./«-i,« iV) ~ iV ~ 1)] 



(A33) 



R{1,N) 



p,(l)p-(2)...p-(jV-l) 
P+(1)P+(2)...P+(A^-1) 



1 + V n™ 



(A34) 



in terms of the ratios (Eq IA30p 



P-ik) 



f, 



k g/3(Jfe + l,lo + 2 — Jfc-l,fc) 



Taking into account that absent links correspond to vanishing coupling Jo.i = = Jn,n+i, one obtains 



p_(lK(2)...p_(iV-l) 



p+{l)p+{2)...p+iN~-l) 



In 



(A35) 



(A36) 
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and 



K[l,I\) -IH — - H T-r — — - 

P+(l)P+(2) 



p_(l)p_(2)...p_(jV-l) 
p+(l)p+(2)...p+(iV-l) 



_|_ h_^PJ2,3 ^ /l_g-0Jl,2+/9J2,3+/9J3,4 _|_ A ^-ff Jl .2 +/3 J3.4 + /3 J4.5 



fl 



-/3Jl,2+/3J4,5+0J5,6 



/3 

A ^-I3Ji.2+0Jn~.2.n^i + 0Jn~i. 



I3Ji,2+I3Jn-i,n 



In-i 

0Jl.2 g0Jl,2+l3J2,a £0^2,3+^^3,4 g0./3 ,4 +/3 J4,; 

— 0Jl.2 



fl 



f'2 



fs 



fi 



/at-I 



N 



= fre-''-^ E 



fu 



fe=i ■'fc 
so that the energy of Eq. IA34I reads 



i?(l,7V) 



N 

E 



fc=l 

i.e. the equihbrium time reads using Eg IA24I 



1 



E' 



^ansatz l^]^-^ 4 j- 

in agreement with the RG resuh of Eq. 11051 derived in the text. 



IN 



(A37) 



(ASS) 



(A39) 



3. Exact renormalization rules in configuration space for escape probabilities Q±{C) 



As a final remark, let us mention the link with previous works concerning renormalization rules in configuration 
space. As explained in [35| . backward master equations satisfy exact renormalization rules in configuration space. 
Upon the elimination of the configuration Co, the surviving configurations C satisfy the same equation as before (Eq 



W,UC)Q+iC) = E ^""(^ ^ C')Q+{C') 



(A40) 



where the renormalized transitions rates evolve with the RG equations 



w. 



Rnew 
out 



(Co) 



(A41) 



These RG rules for backward master equations are exact and can be used '35] , but only for small sizes as a consequence 
of the exponential growth of configurations. The RG rules of Eq. I A4 II have been first derived via a Strong Disorder 
RG approach [36]. 

Note that the RG rules of Eq. IHT] can be rewritten directly for the renormalized probabilities 

.«(C.,.i!i^i^ (A42, 



that evolve according to 



Rne^rr^r, _ W^'-"'(C ^ C) ^^{C) + Trg (Co)7r« (C) 



l-7r«(CoX(C) 



(A43) 
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Appendix B: Dependence on the choice of the dynamics 

In this Appendix, we describe how the equihbrium time depends on the choice of the dynamics satisfying detailed 
balance 

1. Case of the Glauber dynamics 

For the Glauber dynamics, one expects that the dynamical barrier coincides with the maximal energy cost on the 
optimal path between the two ground states. For instance for the one-dimensional random ferromagnetic chain, the 
result of Eq. 11061 satisfies 

1 1, [,o...e,.(^)] ^ (2J,,,0 ^ ^max^ (ur-'^ - C/-) (Bl) 

where U^^'^ '^^ represents the energy of the configuration where the first k spins are (—1), whereas all others spins 
are (+1). 

To better understand the differences with the simple dynamics described below, it is useful to write the result of 
Eq. 11061 for the two smallest sizes, with iV = 2 and = 3 spins 

^glauber (^V = 2) = ^ [l + e^'^'^i'^^] 

^glauber (TV = 3) = ^ [l + 6^'^'^''' + e'^'^'''-^] (B2) 

2. Case of the simple dynamics 



For the 'simple' dynamics, the correspondence of Eq. IBll between the dynamical barrier and the maximal energy 
cost of a single domain wall does not hold, as can be seen already for the one-dimensional case with N — 2 and N = 3 
spins since Eq [94] reads 

t^™P'"(iV = 2) ^^ef^^'-^ 

t^;'"P'"(iV = 3) = I e'^-'^-' + e'^'--^^-'+-^^-''> + e/^-^^-^ (B3) 

For N = 2, the difference by a factor of 2 between the dynamical barriers can be understood from the differences 
between the transitions rates for the simple dynamics (Eq I24p 

simple ^ _ yysimp;e(_|_^ ^ _ g-/3Ji,2 

\YSimple^_^ ^ ^ yysimp;e(_| ^ ^ ^+l3Jl,2 (34-) 

and for the Glauber dynamics feq [^5]) 

p-/3J"i,2 p-2/3Ji,2 



e 



-PJl.2 I + g-2/3Ji 



y^^Glauber^^_ ^ ^ ^Gia«ber(_^_ ^ ) ^ ^ (35) 

For N ^ 2 spins, the equilibrium time is determined by the rate — > H — ) to create a domain-wall when 
starting from one ground state (the time to eliminate the domain- wall is then negligible), and these two rates are 
respectively of order e~^'^'--^ for the simple dynamics and of order e^^^'^^-^ for the Glauber dynamics. One could 
argue that the Glauber dynamics is more 'physical', in the sense that all transitions rates remain bounded near zero- 
temperature, whereas in the 'simple' dynamics transition rates corresponding to a decrease of the energy diverge near 
zero temperature. Nevertheless, one expects on physical grounds that the difference between the dynamical barriers 
of the two dynamics remains of order 0(1), as found in this article for the one-dimensional case and for the tree case, 
and as found in the companion paper |34| for the Dyson hierarchical Ising model. 
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